1 


10/18/10 


2 

3 

4 

5 

6 

7 Global Energy and Water Budgets in MERRA 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


Michael G. Bosilovich * NASA/GSFC/GMAO 
Franklin R. Robertson NASA/MSFC 
Junye Chen U. Maryland/ESSIC and GMAO 


NASA/Global Modeling and Assimilation Office 
Code 610.1 Goddard Space Flight Center 
Greenbelt MD 

Submitted to MERRA Special Issue of Journal of Climate 


Corresponding address : Global Modeling and Assimilation Office, Code 610.1 NASA Goddard Space Flight 
Center, Greenbelt, MD 20771, Michael.Bosilovich@nasa.gov 



1 


1 ABSTRACT 

2 Reanalyses, retrospectively analyzing observations over climatological time scales, 

3 represent a merger between satellite observations and models to provide globally continuous data 

4 and have improved over several generations. Balancing the Earth’s global water and energy 

5 budgets has been a focus of research for more than two decades. Models tend to their own 

6 climate while remotely sensed observations have had varying degrees of uncertainty. This study 

7 evaluates the latest NASA reanalysis, called the Modem Era Retrospective-analysis for Research 

8 and Applications (MERRA), from a global water and energy cycles perspective. 

9 MERRA was configured to provide complete budgets in its output diagnostics, including 

10 the Incremental Analysis Update (IAU), the term that represents the observations influence on 

11 the analyzed states, alongside the physical flux terms. Precipitation in reanalyses is typically 

12 sensitive to the observational analysis. For MERRA, the global mean precipitation bias and 

13 spatial variability are more comparable to merged satellite observations (GPCP and CMAP) than 

14 previous generations of reanalyses. Ocean evaporation also has a much lower value which is 

15 comparable to observed data sets. The global energy budget shows that MERRA cloud effects 

1 6 may be generally weak, leading to excess shortwave radiation reaching the ocean surface. 

17 Evaluating the MERRA time series of budget terms, a significant change occurs, which 

18 does not appear to be represented in observations. In 1999, the global analysis increments of 

19 water vapor changes sign from negative to positive, and primarily lead to more oceanic 

20 precipitation. This change is coincident with the beginning of AMSU radiance assimilation. 

21 Previous and current reanalyses all exhibit some sensitivity to perturbations in the observation 

22 record, and this remains a significant research topic for reanalysis development. The effect of the 

23 changing observing system is evaluated for MERRA water and energy budget terms. 
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1 1. Introduction 

2 In the study of the Earth’s climate, quantifying global water and energy cycling rates and 

3 the associated physical processes more accurately is critical to understanding the climate and its 

4 mechanisms of variability and change from global to local scales. The sun heats the atmosphere 

5 and the surface, thus driving many processes including the transfer of energy and water and 

6 ultimately dynamical transports of these quantities. Trenberth et al. (2009; henceforth TFK09) 

7 provide discussion on the primary water and energy transfer processes, as well as recent 

8 quantitative assessments of various observational data and uncertainties. Even though the TOA 

9 radiative fluxes (please see the appendix for a list of all acronyms) likely have the smallest 

10 uncertainties (order 5 Wm' bias and perhaps an order of magnitude less in precision) refining 

1 1 this observational record is still an active area of research (Loeb et al. 2009). Validating the 

12 observed water cycle observations through global balance shows that the uncertainties are a 

13 fundamental issue (Schlosser and Houser, 2007). While some observational uncertainty is 

14 steadily narrowing, few of the processes have adequate observational representation. In this case, 

15 modeled estimates of the energetics derived from retrospective-analyses (or reanalyses) have 

16 been used to fill gaps in the data. But of course, models themselves, which represent our 

17 understanding of the Earth’s processes, are limited by computational resources and simplifying 

18 assumptions. Models have their own uncertainty and can evolve their own climate, leading to 

1 9 distinct bias when compared with available observations. Data assimilation can produce analyses 

20 of the observed state that constrains the model’s physical results. Retrospective-analyses (or 

21 reanalyses) apply data assimilation across climate time scales in an effort to provide globally 

22 continuous and consistent climate data that exploit both observations and models. However, 

23 individual reanalyses still must contend with uncertainty, where for example, different reanalyses 
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1 respond to global forcing with different circulation perturbations (Chen et al. 2008a). With 

2 several generations of reanalyses to consider, the various data sets generated from these efforts 

3 show large variance in the processes of the global water and energy budgets (Chen et al. 2008 a 

4 and b, TFK09, Bosilovich et al., 2008 and 2009). 

5 Kalnay et al. (1996), Uppala et al. (2005) and Onogi et al. (2007) provide some of the 

6 most important overviews of existing long global reanalyses. In this study, we evaluate the 

7 global energy and water cycles of a new reanalysis, NASA’s Modem Era Retrospective-analysis 

8 for Research and Applications (MERRA, Rienecker et al. 2010). MERRA data are derived from 

9 the Goddard Earth Observing System (GEOS) version 5 data assimilation system, which is a 

10 combination of a NASA general circulation model (Rienecker et al. 2007) and the Gridpoint 

1 1 Statistical Interpolation (GSI) analysis developed in collaboration with the National Centers for 

12 Environmental Prediction (NCEP). During the validation of GEOS-5 and preparations for 

13 MERRA, special attention was given to the water and energy cycles; however, the validation 

14 experiments themselves were limited in time (e.g. Bosilovich et al. 2008). Given the established 

1 5 biases in space and time of the water and energy cycles in existing reanalyses, MERRA water 

16 and energy cycles need to be diagnostically characterized in comparison with the available 

17 observations and reanalyses. 

18 

19 2. Data 

20 a. Reanalyses 

21 There exist several atmospheric reanalyses for the period of 1979 through 2005. The 

22 Japanese 25 year Reanalysis (JRA-25) is quite recent, released for use in March 2006 (Onogi et 

23 al., 2005, 2007), the 45 year ECMWF reanalysis (ERA-40, Uppala et al. 2005), which stops in 



4 


1 August 2002 and the National Centers for Environmental Prediction (NCEP) second reanalysis 

2 (NCEP-DOE, Kanamitsu et al. 2002) represent the second generation of reanalyses. More 

3 recently, ECMWF has also released a short period (1989-present) Interim reanalysis with their 

4 latest data assimilation system (Simmons et al., 2007; Uppala et al., 2008). The NCEP CFSR 

5 (Saha et al. 2010) became available in early 2010. While all these reanalyses assimilate 

6 observations over the recent climate record, new reanalyses will continue to be produced because 

7 of updated model physical processes, enhanced data assimilation methods, increased availability 

8 of computing resource, growing types of observations available for assimilation and improved 

9 observational quality control. 

10 While the data assimilation and numerical model components of the reanalysis system are 

11 fixed for the processing of the climate period (as proposed by Bengtsson and Shukla, 1988; 

12 Trenberth and Olson, 1988), the observing system changes greatly in time. A major change to the 

13 observational record concerns the onset of routine operational remote sensing soundings, starting 

14 with TIROS-N followed by NOAA-6 in 1979, where primarily conventional (e.g. radiosonde and 

15 surface stations) were previously available. Bengtsson et al. (2004) studied the full 40 year time 

16 series of ERA-40 and found that unrealistic trends occurred related to increased satellite 

17 observations in 1979. However, even during the modem satellite observing period, new satellites 

1 8 and measurements start and end at irregular intervals, and ultimately, any given satellite has an 

19 expected lifetime on the order of 10 years, much shorter than what is needed for climate studies. 

20 So, there are numerous changes to the remotely sensed data record (Rienecker et al. 2010) that 

21 involve not only calibration (bias) but also measurement sensitivity and sampling density. These 

22 variations during the satellite era can lead to systematic changes in the reanalysis time series. For 
example, the JRA-25 precipitation record is sensitive to the availability of microwave total 
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1 column water retrievals (Onogi et al. 2005; Bosilovich et al. 2008) beginning in 1987 with the 

2 SSM/I operational satellites. 

3 While analysis state variables are most closely related to observations, the variability of 

4 the physical processes and fluxes among reanalyses can be substantial (Bosilovich et al 2009). 

5 Much of the reanalysis data that provide information about the Earth’s water and energy budgets 

6 come from the model physics, which has been categorized as being closely related to the 

7 numerical model as opposed to the analyzed state fields (Kalnay et al. 1996). Despite the 

8 shortcomings, a particular advantage of reanalyses in climate studies is the availability of all or 

9 many of the Earth’s energy and water budget component fluxes. For example, NCEP reanalyses 

10 have played a significant role in the development of global merged observational precipitation 

11 (CMAP, Xie and Arkin, 1996) and ocean evaporation (Yu and Weller, 2007). Ultimately, 

12 activities such as hydrologic applications would like to make use of reanalyses, but the accuracy 

13 of model physics that control the fluxes requires further development (e.g. Maurer et al. 2001). 

14 The physical terms of the reanalysis budgets generally do not balance even over long periods, 

15 because the atmospheric data assimilation provides additional constraint (or forcing) in the 

1 6 balance of the output data. This is ultimately presented as a residual term in many studies (Roads 

17 and Betts, 2000; Roads et al. 2002). This term, referred to here as the analysis increment, reflects 

1 8 the observations affect on the analysis, and so, as the observations change so do the forcing and 

1 9 the physical response of the model to the forcing. With this term quantified, the budgets can be 

20 studied closely, and some work has used the information to apply corrections to the physical 

21 terms (Schubert and Chang, 1996; Bosilovich and Schubert, 2001; Robertson et al. 2010). 

22 In reviewing the observed global energy budget, TFK09 also compared the reanalysis 

23 energy budgets (specifically, ERA-40 NCEP-DoE R2 and JRA-25), and some similar biases are 
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1 evident. Firstly, the net TOA energy did not balance well, with too much upward flux. However, 

2 JRA-25 bias is related to too much OLR, while NCEP-DoE R2 is due to too much reflected 

3 shortwave radiation, but both imbalances were on the order of 10W/m' . Also, all reanalyses had 

4 excessive evaporation and precipitation leading to stronger global hydrologic cycling. One aspect 

5 of the reanalyses budgets not addressed by TFK09 is the atmospheric imbalance related to the 

6 analysis increment. The observations can act as a source or sink of water and energy, and in their 

7 study, the assimilated observational analysis generally add energy to the system though it is 

8 dissipated in different ways. 

9 b. MERRA 

10 MERRA is a reanalysis of the satellite era (1979-present) using the GEOS-5 data 

1 1 assimilation system. Rienecker et al. (2007) thoroughly describe the MERRA/GEOS-5 numerical 

12 model and data assimilation system, while Rienecker et al (2010) describe the MERRA project. 

13 In addition to the conventional observations (radiosonde, station, aircraft, ship), SSM/I radiances 

14 and various sensor’s retrieved winds, TIROS Operational Vertical Sounder (TOVS) radiances, 

15 Atmospheric Infrared Sounder (AIRS) radiances, and scatterometer wind retrievals (to name a 

16 few) are also assimilated. The experiments were initialized with spun up states from a long 

17 climate model simulation, then 2 years of coarse real data assimilation followed by at least 4 

18 years of real data assimilation (i.e. using the actual observations at the native horizontal 

1 9 resolution, not spinning the same year over again) at the MERRA native grid (Rienecker et al. 

20 2010 describes the spin up experiments). The native MERRA grid has a spatial resolution of V2 0 

21 latitude by 2 A° longitude with 72 hybrid model levels in the vertical (Suarez et al. 2010). The 

22 analysis is performed by the NCEP Gridpoint Statistical Interpolation (GSI; Wu et al 2002). The 

23 model is then updated with an additional model segment that includes an Incremental Analysis 



7 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 


Update, IAU, in the budget equations (Bloom et al. 1996). The shock of the analysis at the 
forecast initialization is greatly reduced, so that the spin down of physical fields (e.g. 
precipitation) is a small factor in this system. This is also where observations affect the model’s 
governing equations (discussed in the next section). Bosilovich et al. (2008) evaluate 
precipitation from the GEOS-5 data assimilation validation experiments for the months of 
January and July 2004 and compare with the existing reanalyses and merged satellite 
observations of precipitation GPCP and CMAP to assess the character of the monthly 
precipitation prior to the production of MERRA. The results were promising but limited due to 
the short period of the validation experiments, and will be briefly revisited later in this study. 
c. Budget Equations 

MERRA output diagnostics encompass all the variables required to balance energy and 
mass budgets for the atmosphere and land. The MERRA enthalpy budget is produced, where 


H = c p T v ( c p is the heat capacity of dry air at constant pressure and T v is the virtual temperature). 
The enthalpy tendency terms are: 

raui rnui rnui 
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Overbars denote a vertical integration over the mass of the atmosphere. The first two terms on 
the right hand side of the equation represent the convergence and release of potential energy. The 
other tendencies represent radiation (RAD), moist processes (MST), turbulent diffusion (TRB, in 
this example, the sensible heat flux). The analysis increment (ANA) is the tendency that is added 
to the prognostic budgets due to the observational analysis. 72 is a small value that in this case 
includes the conversions of energy to/from Kinetic energy (in diffusion and mechanical 
generation), gravity wave drag, and a residual that results from maintaining energy balance in the 
presence of numerical dissipation, each of which are also included in the output diagnostics 
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1 (Suarez et al., 2010). The vertically integrated radiation term can be expanded to its top of the 

2 atmosphere and surface boundary conditions: 

3 \¥\ =(SW t -SW s )-{OLR + LW s ) (2) 

L Ol JRAD 

4 The radiation term includes solar (net shortwave, SW, at the top of the atmosphere, T, 

5 and the surface, S) radiation, the net surface longwave radiation (LWs), and the outgoing LW 

6 radiation (OLR) at the top of the atmosphere. The MST term includes all the heating due to moist 

7 processes, including the condensation heating and evaporation in all phases, and here, the 

8 vertical integral of the MST term is the latent heat resulting from the production of precipitation. 

9 The MERRA vertically integrated total water vapor (w) budget for all phases can be 
10 written as: 

u ^ = -V-(^0 + (£-P) + [f]^+F (3) 

12 w = q v + qi + qi (4) 

13 The change of total water is related to the dynamical convergence of water, and the physical 

14 processes of evaporation and precipitation (sum of convective, large-scale and frozen forms). In 

1 5 the MERRA system, two non-physical terms affect the moisture budget. F represents a very 

1 6 small amount of negative filling, ensuring positive water vapor content. However, the ANA term 

17 represents the analysis increment of water vapor, which is on the order of magnitude of E-P. As 

1 8 MERRA system cycles in time, it performs a forecast, analysis, and then assimilation segment. 

19 The assimilation segment is essentially a model forecast that includes the ANA or analysis 

20 increment tendencies discussed previously. These budgets, and most MERRA output 

21 diagnostics, are derived from the assimilation segment. A result of cycling the system in this way 

22 is that the observational analysis tendencies can be quantified along-side the physical model data. 

23 The result is that long term globally averaged water balances not just E and P, but the analysis 
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1 increment must be accounted for as well. The analysis increment is the forcing applied to a cycle 

2 of the model simulation, developed from the analysis of observations and comparing the analysis 

3 to a forecast cycle (Bloom et al. 1996). Additional information on the formulation of the budgets 

4 is discussed by Rienecker et al (2007) and Suarez et al (2010). 

5 3. Water and Energy Budgets 

6 a. Climatology 

7 TFK09 collect the global energy budget data from various sources, observational and 

8 reanalyses, and close it with consistency arguments from data set intercomparisons, in order to 

9 determine estimates for principal energy flux components and balance. However, each term 

10 exhibits large variations among the different observing systems and reanalyses, so that any 

11 determination of the global average energy budget still includes significant uncertainty. The 

12 spatial and temporal variations then are that much more difficult to know with certainty. We will 

13 characterize the MERRA global energy budget in terms of the existing data and analysis, to 

14 identity those aspects which are realistic at global scales, with some further analysis of regional 

15 scales. 

16 Table 1 compares the results of the TFK09 evaluation of energy fluxes with MERRA 

17 values during the same March 2000 to May 2004 period (see TFK09 for the values from other 

18 reanalyses and observed data, their Table 2). At the top of the atmosphere, the MERRA net 

19 radiative flux is slightly negative (-0.2 W m' ), which compares favorably with an adjusted 

20 ISCCP estimate of 0.9 W m' 2 warming (TFK09). We note that the current unadjusted CERES 

21 TOA imbalance is in excess of 6 Win' 2 (TFK09). The MERRA net flux is closer to the observed 

22 estimates than the JRA-25 and NCEP-NCAR reanalyses. Comparing the TOA components of the 
net radiation, MERRA OLR is larger than the TFK09 estimates (and also observations in their 


23 
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1 Figure 2) while the reflected shortwave radiation is underestimated, suggesting the effect of 

2 clouds is weaker than expected. Similarly, at the surface, downward longwave radiation is 

3 underestimated while the surface downward shortwave radiation is overestimated. 

4 Since MERRA and the reanalyses considered in TFK09 use prescribed SST, the ocean 

5 temperatures and heat content do not respond to the net downward flux at the ocean surface (for 

6 example, 13.8 Wm' 2 in MERRA). While MERRA has a significant global average flux of heat 

7 from the atmosphere to the ocean, the JRA surface flux sign is reversed (TFK09 Table 2). The 

8 NCEP reanalysis shows little average net flux (TFK09, Table 2), especially over the oceans, but 

9 at the expense of unrealistic solar reflection. A significant component of this difference between 

10 MERRA and JRA is the ocean surface evaporation, where MERRA is much lower than JRA. 

1 1 TFK09 (their Table 2) provides merged observation data from HO APS and WHOI evaporation, 

12 but a large discrepancy exists between them. This emphasizes the significant remaining 

13 variations among reanalyses and recent satellite-based estimates. It is anticipated that the 

14 emerging Seaflux data set (Clayson et al., 2010) will narrow the uncertainties for observed data 

15 products. A similar conclusion can also be reached for land fluxes (Vinukollu et al. 2010), where 

1 6 analyses fluxes are collectively biased and more variable than observed products. 

17 

18 b. Spatial variations 

19 As discussed in the previous section, biases in cloud radiative effects plays a role in the 

20 MERRA representation of the surface energy budget. Figure 1 compares the TO A longwave 

21 cloud effect of MERRA to that computed from the SRB data set (Stackhouse et al. 2010). 

22 Positive values indicate that MERRA cloud effect is weaker than SRB. So that, in MERRA, the 

23 tropics have more cloud effect similar to the inter- America seas region (this tends toward boreal 
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1 summer, not shown). On the other hand, extratropical regions have weaker cloud effect 

2 (generally during the boreal winter, not shown). The center of the ITCZ in the eastern tropical 

3 Pacific Ocean is also weaker in MERRA than SRB. The NCEP reanalyses show reasonable 

4 comparison to the SRB data, except that the CFSR SPCZ has a notably weak bias. The ECMWF 

5 and JRA reanalyses have generally weaker cloud effect everywhere. At this scale, high latitude 

6 clouds and radiative effects are not apparent, but Cullather and Bosilovich (2010 a and b) 

7 examine the high latitude water and energy budgets. 

8 A distinctive feature of the MERRA global comparisons in the previous section is the 

9 surface evaporation, especially over the ocean. Where most reanalyses show higher than merged 

10 observation oceanic evaporation, MERRA is lower. Figure 2 shows the surface evaporation 

11 comparison to the WHOI OAFLUX data set (Yu and Weller, 2007). The MERRA low bias is 

12 generally related to extratropical evaporation near the warm western boundary current regions. 

13 The MERRA field has a very strong gradient around the evaporation related to the warm 

14 currents, and the feature is generally smaller scale than in OAFLX (not shown) so that the 

1 5 underestimate in the maximum evaporation near the warm currents may be more related to the 

1 6 differing spatial scales of the data. Away from the warm currents and near the high biases in the 

17 tropical eastern Pacific Ocean, MERRA appears just generally lower than other reanalyses, 

1 8 which consistently overestimate ocean evaporation to varying degrees. 

19 Figure 3 shows the seasonal precipitation comparison among MERRA and other 

20 reanalyses with GPCP (version 2.0 Adler et al., 2003) merged precipitation data. For MERRA, 

21 the tropics precipitation tends to be overestimated, with the midlatitudes underestimated, which 

22 is typical among reanalyses. MERRA’s tropical bias, though, is less than the other reanalyses. 
Many of the precipitation biases apparent here, are also consistent with the biases noted in 
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1 comparing MERRA and TOA longwave cloud effect (Figure 1). For example, in the IAS and 

2 tropical western Pacific regions where the MERRA cloud effect exceeds SRB, the precipitation 

3 is overestimated. Likewise in the southern midlatitudes, MERRA precipitation is slightly 

4 underestimated where the cloud effect is weaker. Reanalyses have been shown to be internally 

5 consistent in terms of cloud-radiation anomalies (i.e. low cloud leads to high OLR), while 

6 differences among the various reanalyses can be substantial (Newman et al., 2000). Several 

7 continental regions’ precipitation are underestimated compared to GPCP, especially South 

8 America. 

9 While MERRA mean precipitation biases are reasonable when compared to GPCP and 

10 other reanalyses (as in Table 1), Figure 4 shows the spatial variability and ultimately the skill of 

1 1 several global reanalyses to reproduce the annual mean distribution of precipitation relative to 

12 GPCP in Taylor diagrams (Taylor 2001; Bosilovich et al. 2008). On these charts, CMAP 

13 provides a secondary reference and sense of the observational uncertainty. Distance from the 1, 1 

14 point represents skill relative to the reference data set (GPCP in this comparison). In the global 

1 5 comparison, MERRA and ERA Interim are the closest to GPCP and CMAP, but both also tend to 

16 be more clustered together (smaller interannual variations of the statistics). CFSR also stands out 

17 ahead of the former generation, though the variance tends to be higher than the current 

18 generation of reanalyses. The global quality in MERRA and Interim can be attributed to the 

19 improvements over ocean regions, but especially the tropical oceans. Figure 5 shows the time 

20 series of the correlation and standard deviation values used to produce the Taylor diagram for the 

21 Globe and Tropics (1 5S-1 5N), as well as the mean bias between the reanalyses and GPCP. While 

22 MERRA has the lowest biases in the tropics and in the group of low bias for the globe, it has an 
increasing trend relative to GPCP. On the other hand, ERA Interim’s high bias in precipitation 
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1 tends to decrease in time. MERRA tropical spatial correlations are higher than any older 

2 reanalysis and comparable to ERA Interim and CFSR. JRA shows a strong change in 

3 precipitation correlation as SSM/I data becomes increasingly available (see also Onogi et al. 

4 2005). The MERRA standard deviation indicates that the variance across the tropics is 

5 comparable to GPCP (which can also be said for NCEP Reanalysis 1 and Interim). The 

6 representation of tropical precipitation in MERRA is much more comparable to GPCP than any 

7 the older reanalyses. At these large scales, the similarity in statistical comparison of MERRA and 

8 Interim with GPCP is remarkable. In the next section, we will evaluate the MERRA time series 

9 of the energy and water budget terms. 

10 4. Time Series 

11 a. Interannual Variability 

12 Previous research has shown that low frequency variability and trends in reanalyses are 

13 problematic and should be considered very carefully. While some fields in certain regions may 

14 be useful indicators of trends (Kalnay and Cai, 2003), there are numerous examples of artifacts 

15 dominating real physical trends from reanalysis data. For example, ERA-40 precipitation has 

16 strong decadal trends and interannual variability of tropical precipitation that do not appear in 

17 observations (Uppala et al. 2005; Andersson et al. 2005). Likewise, JRA-25 precipitation 

18 exhibits a stepwise shift when SSM/I retrieved total column water becomes available for 

19 assimilation (Onogi et al. 2005). In this section, we are evaluating the full time series of MERRA 

20 as an extension to the 5 -year global climatological averages of the energy budget (in the previous 

21 section) and to characterize the system as it changes in time. 

22 The increasing trend in MERRA precipitation bias noted in Figure 5 is separated into 

23 global, land and ocean components in Figure 6a. Precipitation over continental regions shows 
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1 some periods of increase, but the global precipitation time series trends correlate more to the 

2 precipitation over ocean. However, the oceanic trend is not linear, and undergoes a strong 

3 transition period in the late 1990s. Globally and annually integrating the water budget (equation 

4 3), the precipitation is balanced both by evaporation and the analysis increment of water vapor. 

5 Some time variations apparent in global precipitation have similarities in both evaporation and, 

6 especially, the analysis increment (Figure 6b). For example, through the early portion of the time 

7 series, evaporation is steadily increasing and exceeds precipitation by roughly 0. 1 mm d' 1 so that 

8 the moisture increment is negative by a roughly similar amount. However, 1999 and 2001 each 

9 show stepwise changes in the analysis increments, which are generally reflected in the 

10 precipitation. Evaporation also responds to this change by decreasing in 1999, but responds only 

1 1 very weakly in 2001. After 1999, the increments act as a moisture source when the new AMSU 

12 instruments on NOAA-15 (Nov 1998) and NOAA-16 (Jan 2001) become available. Robertson et 

13 al. (2010) discuss the role of the AMSU-A window channels on water vapor and the possible 

14 bias correction uncertainties as the basis for the change in the water vapor increment. 

1 5 Precipitation is clearly sensitive to the changing observing system since the moisture increment 

1 6 changes propagate through the GEOS-5 moisture and radiative physics. 

17 Figure 7 examines the extent of the effect of the satellite instrument changes on the 

18 oceanic surface energy budget. The downwelling shortwave radiation (Figure 7) follows the 

19 water vapor increments (Figure 6b) most closely. In general, the surface shortwave radiation is 

20 decreasing in time while the downward surface longwave radiation is increasing, reflecting the 

21 increasing water vapor and clouds (not shown). The net heating of the ocean surface decreases 

22 substantially over the period with the smallest net heating occurring during the most modem part 

23 of the satellite era (the lowest point occurs in 2001 where it increases steadily after then, Figure 
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1 7). The downward trend in net flux to the ocean surface is almost 10 Wm' over the 30-year 

2 period and if not accounted for could cause significant problems if used for ocean model forcing. 

3 While downwelling shortwave radiation shows a distinct series of changes much like the water 

4 vapor increments or precipitation, the net ocean flux varies more slowly without the sharp 

5 transitions, owing to the combined influences of the surface fluxes. 

6 Figure 8 separates the net radiative flux contributions by LW and SW components at the 

7 surface and top of the atmosphere. The surface net shortwave radiation variations exhibit the 

8 jumps in the times series much like the analysis increments, with the surface longwave radiation 

9 changing in the opposite direction and with smaller amplitude. Likewise, the TOA net shortwave 

10 radiation has much the same variability as surface; both change from positive to negative when 

1 1 AMSU is being assimilated (Figure 8) as cloud effects accompanying precipitation increases 

12 reflection. In the global average, model simulated precipitation heating would be nearly balanced 

13 between the net radiative flux divergence between the TOA and surface. However, for the 

14 reanalysis in the presence of data assimilation, the heating increments must also be considered. 

15 Figure 9 shows the analysis increment of heat in the atmosphere (dHdtANA). This forcing in the 

1 6 global energy budget is ~1 8 W m' in the early part of the period, when satellite observations are 

17 less abundant. After 1998, the heating increments are much lower, ~6 W m' , but clearly 

1 8 sensitive to the availability of AMSU. The other striking point of Figure 9 is the strong negative 

1 9 correlation between the heating increment, and the vertically integrated latent heating (and so the 

20 precipitation) in the atmosphere. Figure 6b shows the water vapor increments (in units of mm 

21 day 1 ), so that there is a strong relationship between the global mean heating increments, water 
vapor increments and precipitation. The magnitude of the response to AMSU in atmospheric 
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2 

1 latent heating and heating increment is ~8 W m' , while the net radiation (surface or TOA) is ~2- 

2 3 W m" 2 . 

3 Figure 6 indicates that the direct effects of the new observations are more consistently 

4 present over the oceans. So that, the possibility exists that these sharp changes could affect the 

5 continental regions through indirect effects, such as moisture transport. Fasullo and Trenberth 

6 (2008a and b) and Trenberth et al. (2009) partition the water and energy budgets (respectively) 

7 into land and ocean components to compute the transport from land to ocean. Table 2 provides 

8 this partitioning for the MERRA water and energy budgets, and solving for the transport. A 

9 fundamental difference from the previously cited calculations is the presence of the analysis 

10 increments in the MERRA budgets. Note that the net radiative fluxes are residuals of very large 

1 1 component flux values. The analysis increments are easily as large as the net radiative fluxes or 

12 the P-E. Considering that the increments represent a corrective tendency accounting for much of 

13 the mean error of the water state at any given time, the magnitude of these terms in the MERRA 

14 budgets is quite large and they cannot be neglected. 

15 For energy, the top of the atmosphere and land surface fluxes appear comparable to the 

16 Fasullo and Trenberth (2008a) estimates. The large downward surface flux over the ocean 

17 mentioned previously is quite apparent in the energy budget, and is likely influencing the ocean- 

1 8 land energy transport (the sign of the flux is opposite to that computed by Fasullo and Trenberth, 

19 2008a). We noted earlier that the imbalance of energy in the JRA-25 ocean radiation is about the 

20 same magnitude as MERRA, but opposite in sign so that different reanalyses could manifest 

21 different water and energy transport characteristics. For the global water budget, the moisture 

22 transport value is generally comparable to that diagnosed by Trenberth et al (2007), and the 

23 magnitudes of the increments are smaller than even the long term average of P-E. This seems 
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1 reasonable for a long term average, but the temporal variations during the period also need to be 

2 considered. 

3 Given the significant extent of the observing system impact on oceanic radiative and 

4 water fluxes, we also computed the time variation of the transport of water and energy between 

5 the ocean and continental areas. Figure 10 (a and b) show the anomalies of water budget terms 

6 for ocean and continental areas, including the moisture transport between them. The shift of the 

7 oceanic water vapor increments is comparable to that of the precipitation. The oceanic water 

8 vapor transport becomes more negative after 1999, indicating more water vapor leaving the 

9 oceanic areas, however, the magnitude of that change is much less than that of the water vapor 

10 increments. On the other hand, for the continental water budgets the variations in precipitation 

11 are much more closely related to the water vapor transport than the increments (Figure 10 b). 

12 The latter is the negative of the ocean transport scaled by land fraction. The heating of the 

13 atmosphere due to the additional water analysis increment is a substantial term in the oceanic 

14 area of the atmospheric heat budget (Figure 10c). However, the variability of this the water term 

15 is nearly opposite of the heating analysis increment, so that when the moisture increment 

16 increases dramatically in 1999, the heating increment decreases. Here, the heat transport is the 

17 sum of all heat transports, including enthalpy, kinetic energy and latent. In this area integration 

18 then, the total transport of heat in Figure 10 c and d does not show a transition in 1999, nor do 

19 any of the component energy transports (not shown). This is also evident in the continental 

20 energy budget, where the increments themselves do not change drastically, but rather, are fairly 

21 steady demonstrating interannual variability anti-correlated to the total heat transport (Figure 

22 lOd). For the whole 30-y period then the land energy budget values are close to those given in 
Table 1 . In contrast to TFK09 these results indicate a weak export of energy from land to ocean 
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1 and a persistent net TOA loss. These losses must be offset by the heating increment to maintain 

2 agreement in temperature and moisture with observations. From a global perspective, this 

3 underscores the fact that the effect of the transition to ATOVS is primarily over the ocean, with 

4 intermittent changes over land occurring though transport processes. 

5 b. Spatial Patterns of Satellite Instrument Change 

6 The previous section shows that in the global average sense, that water vapor and heating 

7 increments change noticeably with the availability of new radiance data (Figure 6b and Figure 

8 10a). The onset of AMSU data availability demarcates distinct climate regimes, from an energy 

9 and water balance point of view. This is especially apparent in the water budget where the 

10 globally averaged water vapor analysis increment changes sign. However, this does not provide 

1 1 more specific information about what regions are affected by the change. 

12 In order to get a sense of the large scale affect of this change in the observing system, we 

13 compare short time averages from before and after the start of AMSU. Decadal time averages 

14 from before AMSU (1990-97) are subtracted from time averages after (2000-08). This may 

1 5 include some real variations in the observed record, but the previous analysis suggests that the 

1 6 analysis increment change (observation system) is the globally dominant factor. Figure 1 1 shows 

17 the change of the water budget component tendencies and fluxes. Precipitation increases in 

18 many places especially in the southern hemisphere midlatitudes, and also in the tropics and 

19 South Pacific Convergence Zone (SPCZ). The increases in precipitation generally correspond to 

20 the locations of changes in the analysis increments, for example, in the southern hemisphere 

21 midlatitudes. However, the SPCZ analysis increments seem smaller than the change in SPCZ 

22 precipitation, but there appears to be an increase in moisture convergence there, so that the large- 
scale circulation has also changed significantly, as there is less convergence (more divergence) in 
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1 the eastern tropical pacific. The changes in evaporation are smaller than the other terms (note the 

2 different contour levels), likely owing to the constraints of prescribed sea surface temperature 

3 and assimilated near surface wind observations. However, the evaporation decreases in the 

4 southern hemisphere midlatitudes where the increased moisture increments have likely decreased 

5 the near-surface humidity deficit in the bulk aerodynamic evaporation term. The increased 

6 atmospheric water vapor has increased the cloudiness, reducing shortwave radiation (increasing 

7 longwave radiation) at the surface. 

8 Figure 12 shows the effect of the AMSU instrument change on the atmospheric dry static 

9 energy budget terms. TOA net (downward) radiation is consistently lower in the presence of 

10 AMSU, driven less by OLR than by increasing reflected shortwave radiation due to increased 

1 1 clouds. The direct effect of heating increments on the radiation components is systematic with 

12 TOA net decreases over persistently cloudy eastern ocean basins, high latitude storm tracks, and 

13 much of Africa. The changes in the analysis increments of heat are quite variable over 

14 continental regions, but they appear to be balanced solely by the vertically integrated heat 

1 5 convergence and release of total potential energy. At the surface, energy losses dominate over 

1 6 ocean regions except in the Southern Ocean where latent heat flux has decreased. However, the 

17 atmosphere also reacts by redistributing the energy dynamically by transport (note the difference 

18 in contour intervals in each panel), consistent with the time series data presented in Figure 10. 

1 9 This may also hold over the oceanic basins as well, except that there is imbalance of energy due 

20 to prescribed SSTs, whereas the land model balances energy over the continental surfaces. The 

21 effect of the AMSU generally reaches all of the water and energy components over oceans, 

especially the warm pool and southern hemisphere midlatitudes, and so these regions will have 
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1 features related to the changing observing system in the long time series. In a companion effort, 

2 Robertson et al. (2010) diagnose these satellite sensor effects in more detail. 

3 In the analysis presented thus far, there appears a widespread systematic change with 

4 time of the water cycle components (Figure 1 1) over the oceans, but more subtle apparent effects 

5 over the continents. However, the tropical continental regions do show some variations, and here 

6 we consider the anomalies in the Amazon River basin and central Africa. Figure 1 3 shows the 

7 time series of root zone soil water for the Amazon River basin (the area is similar to that used in 

8 Bosilovich and Chem, 2006, their Figure 1). In conjunction with the start of AMSU 

9 instrumentation, the soil water becomes systematically higher than in previous years, which 

10 agrees with the increase in precipitation noted in Figure 11a. Figure 14 compares the mean 

1 1 annual cycle of precipitation in the Amazon Basin before and after the start of AMSU for both 

12 GPCP merged gauge-satellite rain rate and that of MERRA. MERRA appears to have a low bias 

13 before AMSU, mostly focused on the transition from dry to wet season. But after AMSU, 

14 MERRA produces more precipitation than the observed data. The AMSU period (1999-2006) 

15 has shifted phase compared to the earlier period, and also compared to observations. It 

1 6 interesting to note that GPCP wet season precipitation in the later period is higher than the earlier 

17 period, which may indicate that there is a real increase in the region. Nonetheless, the change in 

1 8 MERRA is more likely the result of the changing observing systems, and its effect reaching the 

19 large tropical river basin. Given that evaporation, moisture convergence and the water vapor 

20 increments are also changing in time (Figure 11), a more thorough analysis of this region, 

21 beyond the scope of the present study is planned. 

22 In central Africa, a significant hydrologic anomaly occurs, and its sign contrasts the 

23 general changes over the oceans and also over the Amazon River basin. In time, precipitation 
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1 and evaporation are decreasing while the area becomes more divergent despite increasing 

2 vertically integrated water vapor increments. Figure 15a shows the comparison of MERRA 

3 precipitation with that of a station near the center of the low anomaly. MERRA precipitation is 

4 fairly comparable to the local observations early in the reanalysis, but exhibits a sharp drop of 

5 local precipitation in late- 1995. Over the region, the area becomes warm and dry, with 

6 precipitation less than observed. 

7 Further investigation of MERRA’s assimilated observations shows that a single 

8 radiosonde station (Bangui, Station ID 64650) is present in the interior of the continent near the 

9 anomaly, while additional stations are only available at the coasts during most of the period. 

10 Figure 15b shows the monthly mean time series of 850 mb specific humidity anomalies from 

1 1 Bangui, compared with MERRA monthly water vapor anomalies (grid point at 4.5°N, 1 8.667°E). 

12 The radiosonde station water vapor drops substantially in time with MERRA precipitation, but 

13 we also see that MERRA’s water vapor analysis tracks the observations. Since this is the only 

14 nearby radiosonde and satellite data are not as prevalent over land, this station’s observations 

1 5 influence the area radiating out for several hundred kilometers. In time, the surface becomes dry 

16 and warm following the unusually persistent low precipitation. Eventually the atmospheric 

17 circulation is affected, and subsidence forms over the region, possibly related to the AMSU 

18 observation, even further limiting the extent of the precipitation. The feature also appears in 

19 operational analyses surface flux fields comparisons with remotely sensed merged flux data 

20 (Vinukollo et al. 2010). Preliminary investigation in the documentation of this station indicates 

21 that the ground station equipment was changed in the mid-1990s (Leopold Haimberger, personal 
communication 16 September 2010). These results emphasize the continuing need for 
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1 consideration of the input observations along side the reanalyses data to best interpret the 

2 reanalysis results. 

3 

4 5. Summary and conclusions 

5 The MERRA representation of the Earth’s water and energy cycle climatology has 

6 several advantages over existing reanalyses but also weaknesses that have affected past 

7 reanalyses as well. The MERRA climatological precipitation field has a small global bias and 

8 improved spatial correlation compared to the existing global reanalyses, especially across the 

9 tropics, and is also very similar to the ECMWF Interim reanalysis in that regard. While 

10 MERRA’s spatial resolution is finer than many of the previous generation reanalyses, the 

11 number of variables produced is extensive. For example, the systems budget equations are 

12 represented entirely, and can be closed. This closure requires the analysis increment term, that 

13 part of the budgets that allows the model to be influenced by observations. The analysis 

14 increments play a crucial role in the present evaluation of the global water and energy cycles. 

15 Characterization of the strengths and weaknesses requires consideration not only of the model 

1 6 physical processes, but, as emphasized here, the nature of the evolving observational record. 

17 Comparing MERRA global climate energy budgets with previous studies and 

1 8 observations (e.g. TFK09), some consistent biases become apparent. Cloud effects and many of 

19 the surface and TOA radiation components suggest that MERRA clouds are optically weaker 

20 than reality (too few or too thin). This allows excessive shortwave radiation at the surface, 

21 especially over oceans. Time series analysis shows a decreasing trend in the shortwave radiation 

22 at the surface, which is a general improvement to the surface energy balance. However, the time 
series also shows that major stepwise changes occur in MERRA, most notably in the oceanic 
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1 precipitation. These changes are clearly tied to the analysis increment, and are concurrent with 

2 the addition of new satellite instrumentation. This is not a new feature in reanalyses, as the JRA- 

3 25 has a clear dependency on the availability of SSM/I, and CFSR and Interim also show 

4 sensitivity precipitation around the time of NOAA 1 5 and 1 6. Nonetheless, this identifies a 

5 general limitation that must be addressed in future reanalyses. 

6 The strongest shifts in MERRA ’s water and energy budgets coincide with the availability 

7 of AMSU instruments in 1999 and 2001. A concurrent study shows strong evidence that the 

8 sensitivity of the system is related to AMSU-A window channels (Robertson et al., 2010); here 

9 we have assessed the broad impact on the water and energy budgets. The direct effect is largely 

10 related to oceanic regions (especially the southern hemisphere mid-latitudes, tropical western 

1 1 Pacific and Indian oceans), though dynamical transports are significant in linking forcing over 

12 ocean to processes over land, but the effect is more pronounced in some regions than others. The 

13 new data leads to more precipitation, total column water, and clouds, with less net radiation at 

14 the surface and less net radiation leaving at the top of the atmosphere. 

15 The changes in the state of the water balance over the ocean differ in important ways 

16 from those over the land. The transport between land and ocean does not correlate with the 

17 analysis increments over land. Rather, precipitation over land is linked most closely to transport 

18 from the ocean. This suggests that satellite epoch changes (e.g. AMSU-A availability) do 

1 9 indirectly affect MERRA water and energy balance over land through altered moisture and heat 

20 transport, where variations over the Amazon River basin are indicative of this response. Even 

21 so, the dominant data type over land, radiosondes, can strongly influence MERRA state 

22 variables, physics quantities and the trends. This can be demonstrated in central Africa, where a 
single station influences a large region’s water and energy cycles. The observational inputs to 
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1 reanalyses are quite sizable, and have improved demonstrably over the years. These iterations 

2 are clarifying the uncertainties that still exist, as well as their origin in both model physics and in 

3 the quality of the input data stream. It is important for individual researchers to be aware of 

4 variations in the observing system in regards to the science objectives of their research. This also 

5 emphasizes the need for continuing research, development and production of reanalyses. 

6 
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16 6. Appendix: Acronyms 

17 AIRS Atmospheric Infrared Sounder 

1 8 AMSU Advanced Microwave Sounding Unit 

1 9 CERES Cloud's and the Earth's Radiant Energy System 

20 CFSR Climate Forecasting System Reanalysis 

21 CMAP NOAA Climate Prediction Center (CPC) Merged Analysis of Precipitation 

22 ECMWF European Centre for Medium Range Weather Forecasts 

23 ERA ECWMF Reanalysis (40 year or Interim) 

24 GEOS-5 Goddard Earth Observing System (Version 5) 

25 GPCC Global Precipitation Climatology Centre 

26 GPCP Global Precipitation Climatology Project 

27 GSI Gridpoint Statistical Interpolation 

28 HOAPS Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data 

29 IAS Inter-America Seas 

30 IAU Incremental Analysis Update 
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1 ISCCP International Satellite Cloud Climatology Project 

2 ITCZ Intertropical Convergence Zone 

3 JRA-25 Japanese 25 year reanalysis 

4 LW Longwave 

5 MERRA Modem Era Retrospective-analysis for Research and Applications 

6 NASA National Aeronautics and Space Administration 

7 NCAR National Centers for Atmospheric Research 

8 NCEP National Center for Environmental Prediction 

9 NOAA National Oceanic and Atmospheric Administration 

1 0 NRA NCEP-NCAR reanalysis 

1 1 OAFLUX Objectively Analyzed air-sea Fluxes 

1 2 OLR Outgoing longwave radiation 

1 3 PW Petawatt 

14 SPCZ South Pacific Convergence Zone 

1 5 SRB Surface Radiation Budget 

16 SSM/I Special Sensor Microwave Imager 

17 S ST Sea Surface Temperature 

1 8 TIROS Television Infrared Observation Satellite 

1 9 TOA Top of atmosphere 

20 TOVS TIROS Operational Vertical Sounder 

21 WHOI Woods Hole Oceanographic Institution 

22 
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1 8. Tables 

2 Table 1 Energy fluxes (March 2000 - May 2004) partitioned by Global, Global Land and Global 

3 Ocean averages and compared to the estimates developed by TFK09. TFK09 also provide 

4 ISSCP-FD, NCEP reanalyses, JRA-25 and WHOI and HOAPS ocean fluxes for 

5 comparison. 


a.) TOA 


Global 

Solar In 

Solar Refl. 

Albedo % 

Solar Net 

OLR 

NET down 

MERRA 

341.3 

99.6 

29.2 

241.7 

242.0 

-0.2 

TFK09 

341.3 

101.9 

29.8 

239.4 

238.5 

0.9 

Land | 

MERRA 

325.7 

106.8 

32.8 

218.9 

232.9 

-14.0 

TFK09 

330.2 

113.4 

34.4 

216.8 

232.4 

-15.6 

Ocean | 

MERRA 

347.9 

96.6 

27.8 

251.2 

245.7 

5.5 

TFK09 

345.4 

97.8 

28.3 

247.7 

240.8 

6.9 


b.) Surface 


Global 

Solar 

Absorb 

Solar Net 

Solar 

reflect 

LH 

SH 

LW Up 

LW Down 

Net LW 

NET down 

MERRA 

72.7 

169.1 

23.6 

76.4 

18.3 

394.3 

330.3 

64.0 

CO 

4*. 

TFK09 

78.2 

161.2 

23.1 

80.0 

17.0 

396.0 

333.0 

63.0 

1.2 I 

Land | 

MERRA 

68.7 

150.2 

46.7 

44.9 

33.2 

364.6 

294.6 

70.0 

2.0 

TFK09 

78.0 

145.1 

39.6 

38.5 

27.0 

383.2 

303.6 

79.6 

0.0 

Ocean ] 

MERRA 

74.3 

176.9 

14.0 

89.5 

12.1 

406.6 

345.2 

61.5 

13.8 

TFK09 

78.2 

167.8 

16.6 

97.2 

12.0 

400.8 

343.2 

57.6 

1.3 
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1 

2 Table 2 The land and ocean water and energy budgets including transport (dwdtDYN and 

3 dHdtDYN) between land/ocean from MERRA 30 years average. The water budget terms 

4 are evaporation (E), precipitation (P), water vapor analysis increment (dwdtANA). The 

5 energy budget terms are net fluxes at the surface (SFC) and top of the atmosphere (TOA), 

6 analysis increment of dry static heat (dHdtANA) and heating of due to the analysis of 

7 water (dHdtANAw) 


Water Budget (xlO 3 Km 3 yr' 1 ) 



E 

P 

dwdtANA 

dwdtDYN 

Land 

87 

116 

-9 

39 

Ocean 

405 

368 

1 

-39 

Energy Budget (PW) 



SFC 

TOA 

dHdtANA 

dHdtANAw 

dHdtDYN 

Land 

-0.1 

-2.0 

3.4 

- 0.7 

-0.6 

Ocean 

-6.5 

2.9 

2.8 

0.1 

0.6 
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9 

10 9. List of Figures 

1 1 Figure 1 Annual differences (1990-2002) of MERRA and other reanalyses longwave cloud effect 

12 from that of the Surface Radiation Budget (SRB) data. Here, the longwave effect is the 

1 3 TOA all-sky outgoing longwave radiation (OLR) minus the clear-sky outgoing longwave 

14 radiation, so that positive difference indicates the reanalysis cloud effect is less than SRB. 

1 5 There global area average of each map is included in the upper right comer of each panel. 

16 Additional figures are available from The MERRA Atlas (GMAO, 2010). 
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Figure 2 Annual differences (1990-2002) between MERRA and other reanalyses ocean surface 
latent heat flux and that of OAFLUX merged data (Weller and Yu, 2007). Positive flux is 
directed upward. 

Figure 3 Annual differences (1990-2002) between MERRA and other reanalyses precipitation 
with GPCP merged data (Adler et al. 2003). 

Figure 4 Taylor diagrams of annual mean precipitation from reanalyses using GPCP as a 
reference and CMAP as an additional observing reference. Each panel shows the 
statistics for different regions: (a) globe, (b) land, (c) ocean and (d) tropics. The red and 
blue lines show limits of expected high and low correlation as determined by comparing 
GPCP and CMAP observations. See Bosilovich et al (2008) for details. 

Figure 5 Time series of annual mean statistics used in the global (a-c) and tropics (d-f) Taylor 
diagrams (Figure 4 a,d) for (a,d) mean difference form GPCP, (b,e) spatial correlation to 
GPCP and (c,f) standard deviation (here the black line is GPCP standard deviation, while 
the Taylor diagrams are normalized to GPCP standard deviation). Averages consider only 
gridpoints where valid CMAP observation only data exists. 

Figure 6 Global annual averages of (a) precipitation including land-only and ocean only 
averages, and (b) the global annual averages of the water budget including precipitation 
(P), evaporation (E) and analysis increment of water vapor (Qvinc). 

Figure 7 Global ocean-only (a) annual anomalies for surface energy budget terms latent heat flux 
(-LE), sensible heat flux (-FIs), downward shortwave radiation (SWdn) and downward 
longwave radiation (LWdn) and (b) annual mean Net downward flux at the ocean 
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surface. The mean values of the terms are oriented positive down to the surface. Mean 

2 values of the anomalies are included in the legend. Units are W m' . 

3 Figure 8 Global annual anomalies for (a) net radiation terms at the surface, shortwave (SWnet), 

4 longwave (LWnet) and the Net Radiation, (b) surface heat fluxes (positive down) 

5 including the Net heating and (c) radiation terms at the top of the atmosphere (TOA) net 

6 shortwave radiation (SWnet), upward longwave (LWup) and net top of atmosphere 

7 radiation (TOAnet). The means for each term are included in the legends. Units are in W 

8 m' 2 . 

9 Figure 9 Energy terms for the moist processes (MST, essentially latent heating due to 

10 precipitation) and the heating due to analysis increments. Units in WnT . 

11 Figure 10 Separate terms for the integrated water mass and energy budgets over land and ocean 

12 to compute land-ocean transport (a) ocean water budget anomalies, (b) continental budget 

13 (see table 2 for the time means and variable definitions), (c) oceanic energy terms and (d) 

14 land energy terms. 

15 Figure 11 Change (2000-08 minus 1990-97, after AMSU minus before) of the atmospheric 

1 6 hydrology budget terms Precipitation (a), Evaporation (b) Incremental Analysis Update 

17 for water vapor (c) and moisture convergence (d). The change of evaporation is smaller 

1 8 than other terms and is scaled different from other terms. 

19 Figure 12 Change (2000-08 minus 1990-97, after AMSU minus before) of the atmospheric 

20 energy budget terms Net TOA radiation (a) Net Surface Energy Flux (b), Incremental 

21 Analysis Update for dry static heat (c) and heat convergence. 
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1 Figure 13 Time series of MERRA monthly mean root zone soil wetness fraction 

2 (nondimensional), area averaged over the Amazon River basin. 

3 Figure 14 Mean annual cycle of monthly precipitation Area averaged over the Amazon River 

4 basin from MERRA before AMSU data (1979-1998, dashed line) and MERAA when 

5 AMSU data is assimilated (1999-2006), solid line, and also the corresponding GPCP 

6 (Adler et al. 2003) basin averaged precipitation (• for 1979-1998, and A for 1999- 

7 2006). Units are mm day' 1 . 

8 Figure 15 Time series of monthly MERRA data compared to observations at the MERRA point 

9 (4.5°N, 1 8.67°E). (a) GPCC precipitation (dots with thin line; Schneider et al. 2008) and 

10 MERRA precipitation (thick solid line) and (b) Bangui radiosonde observations (dots, 

11 station ID 64650, 4.24°N, 18.31°E) of 850 mb specific humidity anomaly compared to 

12 MERRA. The anomalies were calculated by removing the mean annual cycle. 
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10. Figures 
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Figure 1 Annual differences (1990-2002) of MERRA and other reanalyses longwave cloud effect 


from that of the Surface Radiation Budget (SRB) data. Here, the longwave effect is the TOA all- 
sky outgoing longwave radiation (OLR) minus the clear-sky outgoing longwave radiation, so 
that positive difference indicates the reanalysis cloud effect is less than SRB. 
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Latent Heat Flux Difference 
(Reanalysis-OAFLUX W m 2 ) 
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2 Figure 2 Annual differences (1990-2002) between MERRA and other reanalyses ocean surface 


3 latent heat flux and that of OAFLUX merged data (Weller and Yu, 2007). 
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3 Figure 3 Annual differences (1990-2002) between MERRA and other reanalyses precipitation 


4 with GPCP merged data (Adler et al. 2003). 
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Figure 4 Taylor diagrams of annual mean precipitation from reanalyses using GPCP as a 
reference and CMAP as an additional observing reference. Each panel shows the statistics for 
different regions: (a) globe, (b) land, (c) ocean and (d) tropics. The red and blue lines show limits 
of expected high and low correlation as determined by comparing GPCP and CMAP 
observations. See Bosilovich et al (2008) for details. 
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Global Precip Bias (mm day -1 ) 
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Figure 5 Time series of annual mean statistics used in the global (a-c) and tropics (d-f) Taylor 
diagrams (Figure 4 a,d) for (a,d) mean difference form GPCP, (b,e) spatial correlation to GPCP 
and (c,f) standard deviation (here the black line is GPCP standard deviation, while the Taylor 
diagrams are normalized to GPCP standard deviation). 
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Figure 6 Global annual averages of (a) precipitation including land-only and ocean only 
averages, and (b) the global annual averages of the water budget including precipitation (P), 
evaporation (E) and analysis increment of water vapor (Qvinc). 


Ocean Surface Flux Anomalies (W m 2 ) ^ 



r-- oo oc oo oo 

^ On On O'* On 


Os Os O'* 


— -LE (-87.7) - - -Hs (-11.7) 

S Wdn (1 95 . 8) — LWdn (342.7) 



CTN’-^c^unr-'-c^'— 

r~-ooooooooooc^<^<^cr'(<^ooooo 


Figure 7 Global ocean-only (a) annual anomalies for surface energy budget terms latent heat flux 
(-LE), sensible heat flux (-Hs), downward shortwave radiation (SWdn) and downward longwave 
radiation (LWdn) and (b) annual mean Net downward flux at the ocean surface. 
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Figure 8 Global annual anomalies for (a) net radiation terms at the surface, shortwave (SWnet), 
longwave (LWnet) and the Net Radiation, (b) surface heat fluxes (positive down) including the 
Net heating and (c) radiation terms at the top of the atmosphere (TOA) net shortwave radiation 
(SWnet), upward longwave (LWup) and net top of atmosphere radiation (TOAnet). The means 
for each term are included in the legends. Units are in W m 2 . 
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Moist Process and Analysis Heating (W nr 2 ) 
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Figure 9 Energy terms for the moist processes (MST, essentially latent heating due to 

5 precipitation) and the heating due to analysis increments. Units in Wm' . 
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Figure 10 Separate terms for the integrated water mass and energy budgets over land and ocean 
to compute land-ocean transport (a) ocean water budget anomalies, (b) continental budget (see 
table 2 for the time means and variable definitions), (c) oceanic energy terms and (d) land energy 
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Figure 11 Change (2000-08 minus 1990-97, after AMSU minus before) of the atmospheric 
hydrology budget terms Precipitation (a), Evaporation (b) Incremental Analysis Update for water 
vapor (c) and moisture convergence (d). The change of evaporation is smaller than other terms 
and is scaled different from other terms. 
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Figure 12 Change (2000-08 minus 1990-97, after AMSU minus before) of the atmospheric 
energy budget terms Net TOA radiation (a) Net Surface Energy Flux (b), Incremental Analysis 
Update for diy static heat (c) and heat convergence. 
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Amazon Basin Root Zone Soil Wetness Fraction 



2 Figure 13 Time series of MERRA monthly mean root zone soil wetness fraction 

3 (nondimensional), area averaged over the Amazon River basin. 
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Figure 14 Mean annual cycle of monthly precipitation Area averaged over the Amazon River 
basin from MERRA before AMSU data (1979-1998, dashed line) and MERAA when AMSU 
data is assimilated (1999-2006), solid line, and also the corresponding GPCP (Adler et al. 2003) 

basin averaged precipitation (• for 1979-1998, and A for 1999-2006). Units are mm day' 1 . 
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3 Figure 15 Time series of monthly MERRA data compared to observations at the MERRA point 

4 (4.5°N, 18.67°E). (a) GPCC precipitation (dots with thin line; Schneider et al. 2008) and 

5 MERRA precipitation (thick solid line) and (b) Bangui radiosonde observations (dots, station ID 

6 64650, 4.24°N, 18.31°E) of 850 mb specific humidity anomaly compared to MERRA. The 

7 anomalies were calculated by removing the mean annual cycle. 
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